Dynamical qubit controlling via 
pseudo-supersymmetry of two-level systems 



Boris F. Samsonov and V. V. Shamshutdinova 

Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russia 

Abstract. For a flux qubit considered as a two-level system, for which a 
hidden polynomial pseudo-supersymmetry was previously discovered, we propose 
a special time-dependent external control field. We show that for a qubit placed 
in this field there exists a critical value of tunnel frequency. When the tunnel 
frequency is close enough to its critical value, the external field frequency may 
be tuned in a way to keep the probability to detect a definite direction of the 
current circulating in a Josephson-j unction circuit above 1/2 during a desirable 
time interval. We also show that such a behavior is not much affected by a 
sufficiently small dissipation. 
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1. Introduction 

In the past few years superconducting circuits based on Josephson tunnel junction 
have attracted much attention both from theoretical and experimental viewpoints 
as possible candidates for the implementation of quantum computer (see, e.g. Rcfs. 
[HHIE])- Usually they represent a small Josephson-junction circuit, called a Cooper- 
pair box, which consists of a small superconducting electrode connected to a reservoir 
via a Josephson junction [3J. For a flux qubit the circuit with a very small inductance 
containing three Josephson junctions is described (in appropriate units) by the 
following two-level Hamiltonian [3J: 

H q = -A<j x - e{t)a z . (1) 

Here A is the tunnel frequency, and s(t) is a time-dependent field (bias) which is 
controlled by an externally applied flux. Although in general A is a function of e, 
it varies on the scale of the single junction plasma frequency, which is much above 
the typical energy range at which the qubit is operated [4]. We thus can assume A 
to be constant for the purpose of this paper. Solving the Schrodinger equation with 
Hamiltonian {T]), 

M{t)=H q V, *=(V»i,^) T (2) 

(superscript "T" means the transposition and the dot over a symbol means the 
derivative with respect to time), we obtain the probability of a definite direction of 
the current circulating in the ring, i.e. P^ — |V>i| 2 is the probability of the clockwise 
direction of the current and = \ip2\ 2 is the probability of the opposite current 
direction. 

One of the most important problem in quantum computation is connected with 
the possibility of controlling the state of an array of qubits. Typically the simplest 
two-qubit operations are generated by interplay of the coupling between qubits and 
local fields. Much theoretical attention has been recently paid to studying the 
controllable coupling between qubits of different types (see Ref. [5] and references 
therein). Recently it has been shown [6] that in the simplest and the most important 
from engineering viewpoint case of an "always on" and fixed coupling, a two-qubit 
Hamiltonian may be decoupled and the control problem is, in particular, reduced to 
finding the evolution of a one qubit placed in a time-dependent external control field. 
This observation shows an additional importance of controlling a one qubit state. This 
is the main subject we devote the present paper. 

Usually, the probabilities, as functions of time, show an oscillating behavior (cf. 
famous Rabi oscillations, see e.g. Ref. [7]). But for some specific external fields this 
character may be changed drastically thus showing a possibility to control the qubit 
state. Up to now such a possibility is known to be mostly related to oscillating external 
fields 0. 

Recently [9] it was proposed to consider a two-level Schrodinger equation as a 
Dirac equation with a non-Hcrmitian Hamiltonian where time plays the role of a 
spatial variable. This possibility revealed a hidden pseudo-supersymmetry which may 
be associated with a two-level system [10] and lead to discovering new time-dependent 
external fields where a two-level system admits solutions in terms of elementary 
functions. An advantage of analytic solutions is the possibility of a careful analysis of 
their properties which may reveal unexpected peculiarities [11] . In this paper we apply 
these results to show the possibility of controlling the qubit state with an external 
field of a special configuration. We show that there exists a critical value of the tunnel 
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frequency A. While the tunnel frequency approaches this critical value, the probability 
P^(t) oscillates around a value exceeding 1/2 with a decreasing amplitude and after 
the critical value is reached it becomes a function monotonously increasing up to 
limiting value equal 3/4. Then using the property that this special excitation regime 
is, in fact, a limiting case of a more general oscillating external field, we demonstrate 
that one can control a definite direction of the current in the ring (i.e. the qubit 
state) during a desirable time interval. Finally we also show that such a behavior is 
not much affected by the presence of a reasonably small dissipation featuring open 
quantum systems. We start with reminding the reader main constructions leading 
to polynomial pseudo-supersymmetry in two-level systems which is done in the next 
section. 



2. Polynomial pseudo-supersymmetry of two-level systems 

Similar to the conventional supersymmetry in quantum mechanics (for recent 
developments see a special issue of 2004 J. Phys. A: Math. Gen. 34 (43)), pseudo- 
supersymmetry is based on intertwining and factorization relations. However, in 
this case Hamiltonians are non-Hermitian and a specific automorphism, that defines 
pseudo-adjoint operators, should be involved [T5J[T3]. In particular if A is a linear 
operator and r\ is linear, Hermitian, invertible operator then 

A* = ri- l A + Ti, (3) 

where '+' sign denotes the usual (e.g. Laplace) adjoint operation, by definition is the 
pseudo-adjoint of A with respect to 77. The operation of formal (Laplace) conjugation 
obeys the standard rules {AB)^ — B'A\ (d/dt)^ = —d/dt and corresponds to the 
transposition of a matrix accompanied by the complex conjugation of its elements. 
Operator B is said to be pseudo-Hermitian with respect to 77 if B* = B, i.e. 

B = rf- x B + r}. (4) 

Basic properties of pseudo-Hermitian operators are discussed in detail in [13] . 

Recently |10j it was observed that a polynomial pseudo-supersymmetry may be 
associated with a two-level system interacting with a classical (i.e. not quantized) 
electromagnetic field. The method is based on the possibility to rewrite the 
Schrodinger equation @ that governs the evolution of the system in the form of 
a one-dimensional stationary Dirac equation 

hoV = EV, (5) 

where the time plays the role of a spatial variable and the Dirac Hamiltonian 

ho = io- x ^+V (t), V (t) = ia y e Q (t), E = A (6) 

is non-Hermitian, with the subsequent application of the well-developed intertwining 
operators technique [14J . Function e(t) plays the role of a "potential". 

It is easy to see that ho (JB]) is pseudo-Hermitian with respect to 77 = a x . The next 
ingredient of the method is based on the existence for any real- valued function £o(t) 
such real- valued function £i(i) and operator L (intertwiner) that 

Lh Q = hiL (7) 

where hi = io~ x d/dt + V\{t) with Vi = ia y Ei(t). The pseudo-Hcrmiticity of ho results 
in the following factrorizations 

r^L+r/L = h\ - A 2 , LnL+r] = h\ - A 2 . (8) 
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The constant matrix A = diag(A, —A) in © is called the (matrix) factorization 
constant (A is also called the factorization constant). Formulas {SJ present a 
generalization of the factorization properties of transformation operators that take 
place in the case of Hermitian one-component Hamiltonians [15) . If now we introduce 
matrix operators (in block-matrix forms) 

s). *-(s jl : j ) « 9 > 

then the intertwining ([7]) and factorization (jSj) relations may be rewritten as the 
following set of commutation and anti-commutation relations between these operators: 

Ql 2 = 0, HQ h2 = Q h2 H, QiQ 2 + Q2Q1 = H 2 - diag(A, A) (10) 

which indicate on the simplest quadratic pseudo-superlagebra. Note that the 
subsequent application of this technique leads to a more general polynomial pseudo- 
supersymmetry |10) . On the other hand if we start from Vq with known solutions to 
the Dirac equation (J5J then solutions of the same equation with V\ can be obtained 
by applying L to the previous solutions. In this way new exactly solvable two-level 
potentials are obtained [10] which we use in the next section to demonstrate the 
possibility of the dynamical qubit controlling. 



3. Dynamical qubit controlling 



* (*) = -eo + tt™ ■ (ID 



Consider first the case when the external control field e = £\{t) changes in the following 
way: 

4e 
1 + 4e§i 2 

Parameter £0 gives us the possibility to choose a suitable time scale since after re- 
scaling r = 2eoi, and redefining parameter A, A — 2eo^, we obtain the Schrodinger 
equation with the Hamiltonian 

H = -Sa x - e(r)a z , c(r) = e x (r) = -\ + — (12) 

2 1 + T z 

for which exact analytic solutions are known [10] . Therefore imposing the initial 
condition P* (0) = we can write down an explicit expression for the probability 

p l {T)= {0 2 -l)(6 2 + A) r 2 , {9 2 -l){e 2 -A) 



29 4 1 + T 2 26> 6 (1 



r 



2) 



x [9 2 - 4 - (6 2 - 4 + 2 t 2 ) cos 9t + A9t sin 6t] (13) 

where we have introduced 8 = yl + AS 2 . From here it is clearly seen that pI(t) is an 
oscillating function provided ^ 2 (5 ^ V3/2). For 9 = 2 (5 — v3/2) equation (fl~3| 
yields 

which is a function monotonously increasing from zero at r = till the value 3/4 
for r 3> 1 (solid (black) line in Fig. [T|). Since 5 differs from A only by the scaling 
factor 2sq we will call S tunnel frequency as well. We note a decrease of the oscillation 
amplitude when S approaches its critical value equal v3/2. This is why for S close 
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Figure 1. (Color online) Evolution of P-j* probability at 5 = V3/2 solid, dot- 
dashed and double-dot-dashed (black, red and green) lines, at dephasing (Y v ) and 



relaxation (IV) rates: F r 
double-dot-dashed (green) line; 
blue) lines resp. 
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Figure 2. (Color online) (a, c, d) Evolution of probability at 8 = %/3/2 + 0.1, 
at dephasing (r^,) and relaxation (r r ) rates: Y v = 0.1, Y r = 0.05 dot-dashed 
(red) line and r r = 0.2 double-dot-dashed (green) line; (a) = 0.105, (b) 
= 0.205 and (c) = 0.314. 




12 3^ 01234,5 



Figure 3. (Color online) <5-dependence of time-averaged probabilities for the 
closed system solid (black) lines and with dephasing (r v ) and relaxation (IV) 
rates: Y r = T v = 0.01 dot-dashed (red) lines; (a) Y r = Y v =0.1 and (b) 
r r = Y v = 0.05 double-dot-dashed (green) lines. 



enough to the value y/3/2 the minimal value of the probability P/(r) for r > 2 exceeds 
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1/2 (see dashed (blue) and dotted (purple) lines in Fig. [J). 

The time-averaged probability as shown in Fig. [3^ by solid (black) line exhibits 
a maximum. Its analytic expression 

T} = 26> 5 + 4§ \ (15) 
(1 + 45 2 ) 2 

allows us to get the exact position of the maximum which is 25/32 ss 0.78 at 
5=^fbJV2. 

The result we have just obtained suggests us to consider a more general case 
[HlUnij where function e(r) being periodical depends on three parameters. One of 
them we fix by re-scaling both time and another parameter (frequency oj appearing 
in Eq. ([To]) ) in a way as it has been done above, thus obtaining the field 

2iu 2 

e 2 (r) = -l/2 ; r (16) 

y ' 1 bcos{2ujT + ip) - 1/2 1 ' 

where b 2 = 1/4 — 2 > 0. The remaining parameter ip can also be eliminated by 
shifting the time origin so that we put ip = thus reducing external fields to a one 
parameter (0) family. It is important to note that the previous result (|12|) may be 
obtained from here at if = arctanw — \ arctan ^ in the limit 0-^0. 

According to [10] the analytic expression for P^ij) reads 

zLA 2 

46 2 b [Q - 0(6 + b 2 - S 2 )9 sin(20r) sin(0r)] 
9 2 (b 2 +S 2 ) 2 (2bcos(20T) - 1) 

where 

Q = b(l + 2b)8 2 cos 2 (|0t) sin 2 (0r) 
+ 40 2 (fe - 2S 2 ) cos 2 (0t) sin 2 (±6>t) . 

To show the possibility to control the qubit state by external field (fT6|) we plot 
function -P 2 ^( T ) f° r ^ close to its critical value and for different values of 0. Since 
the above considered case corresponds to = 0, we show in Fig. [2^i (solid (black) 
line) its behavior for = 0.105 which is rather close to zero. During a sufficiently 
long time interval the probability oscillates between 0.6 and 0.8 after which it falls to 
zero. The closer gets to zero, the longer this period becomes and the closer to 1/2 
becomes, the more this interval is reduced (see the solid (black) lines in Figs. [2^,, 
[2jo and[2j:). We have to note that in the limit — > 1/2 (b — > in Eq. (fl6|) 1 function 
C2(t) tends to a constant value equal 1/2. This signal reproduces the Rabi oscillations 
with the frequency 2^JS 2 + 1/4. Thus, may be considered as a continuously tunable 
parameter of the external field (fT6|) with the help of which, starting with the usual 
Rabi oscillations, one may fix the clockwise current direction as long as desirable. We 
would like to stress that the range of 5, for which the probability exceeds the value 
1/2, is rather large, i.e. 0.6 < S < 1.1. This may facilitate its experimental detection. 

Analytic expression for (t) shows us that it represents a complicated 
superposition of two oscillating functions with frequencies 8 and 20. Therefore if these 
frequencies are close enough to each other one may observe a beating phenomenon 
(Fig. [4]). In this case the oscillations with the small amplitude and the frequency 
close to the Rabi frequency take place at the background of the oscillations with the 
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Figure 4. (Color online) Time dependence of clockwise current direction 
probabilities at = 0.49 and <5 = 0.1 solid (black line), at dephasing (r^,) and 
relaxation (IV) rates: T r = Y v = 0.01 dot-dashed (red) line and T r = F v = 0.02 
double-dot-dashcd (green) line. 



amplitude close to 1 and the very small frequency denned by the difference between 
20 and 8. 

Another aspect we would like to emphasize is that this type of the external field is 
not unique. As we show below, there exist other possibilities for the time dependence 
of the external field exhibiting a similar feature. 

Consider an exactly solvable model with a bit more complicated form of function 
e(r) = e 3 (r) [10] 

*{t) = -\ +7?-^ 4 + 6r 2 -3) (17) 

where Q = r 6 + 3r 4 + 27r 2 + 9. The clockwise current direction probability P^(t) = 
P 3 (r) in this case reads 

P3 (r) = 4( ^ 8Q ^ )r2 [144(1 + r 2 ) + <? 4 (r 2 + 9) 2 - 240 2 (5r 2 + 9)] 

+ (g2 g r Q Ql [Q2 sin 2 (±r#) + g 3 sm(r0)] (18) 



where 



Q 1= [(0 + l) 2 -5)] [(0-l) 2 -5)] 

Q 2 = 6» 4 Qo + 144(1 + r 2 ) - 126> 2 (5t 4 + 6r 2 + 9) 

Q 3 = 68t [0 2 (t 4 + 2r 2 + 9) - 12(1 + t 2 )] . 

We note that this is just the second term in the right hand side of Eq. (fl8|) which 
is responsible for time oscillations. Therefore if Q\ = the oscillations in the time 
dependence of the probability disappear and once again it acquires a monotonous 
character. But this time since Q\(0) is a bi-quadratic function, in contrast to the 
previous case, the clockwise current direction probability turns from an oscillating to 
monotonous character at two possible values of parameter 9, 8 = y/E ± 1. In these 
cases the behavior of probabilities P^(t) and p|(t) is illustrated in Figs. [5]a, and 
\5jp (solid (black lines)) respectively. We thus observe for P 3 ^(r) an effect similar to 
that described above for Pj *(r) and, in a sense, the opposite behavior of P3 (r). The 

existence of the critical value is reflected also by the averaged probability P 3 which is 
plotted in Fig. (solid (black) line). It also has a simple analytic expression 

yi _ 2 13-8o 2 + 16^ 
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Figure 5. (Color online) Probabilities (a) Pg(r) at 8 = \/5 — 1 solid (black) 
line, 9 = v6 — 1.1 dotted (violet) line and 8 = v5 — 0.9 dashed (blue) line; (b) 
P 3 T (t) at = v / 5 + 1 solid (black) line, d = a/5 + 1-5 dotted (violet) line and 
= + 0.5 dashed (blue) line. 



with a maximum P 3 « 0.91 at £ w 0.34. 

The next question we study is how the effect, observing for an idealized closed 
system, is influenced by a dissipation featuring open quantum systems [16 . To 
make rough estimations we use a phenomenological approach in the density matrix 
formalism (see e.g. |17]V In this approach a weak coupling of a system to the 
environment is described by two parameters, dephasing T v and relaxation T r rates. 
Under the initial condition 2p(0) = I with / being identity matrix the elements of the 
density matrix (cf. 18J ) 

1 + Z X-iY ' 
X + iY 1 - Z 

satisfy the Bloch equations [T7] 

X = -2e{t)Y-T v X 

Y = -25Z + 2e(t)X - T V Y 

Z = 2SY - r r (Z - Z(0)) . 

Probabilities pJ"T are defined only by the diagonal entries of the density matrix, 
= (l T Z(t))/2. 

The relaxation and dephasing effects are shown in Figs. HHUby dot-dashed (red) 
and double-dot-dashed (green) lines. As expected, they disturb the system. The 
influence of dephasing is more crucial and it should not exceed 5 percent of 5 value for 
bias l|16p. We observe an interesting phenomenon concerning the relaxation. When the 
probability falls to zero the relaxation smoothes this behavior and can even revert it 
(see green line in Fig. [5^). Thus, in the current case the relaxation may be considered 
as helping to keep the state of the qubit unchanged. 

4. Concluding remarks 

One of the most popular way to control the qubit state is to drive a two level 
system with microwave pulses (see e.g. |19j ) where, in general, the probability P^(t) 
displays (possibly controllable) Rabi oscillations. On the other hand it is known that 
with a specific monochromatic driving force one can "freeze" the state of a two-level 
system (so called "coherent destruction of tunnelling" [8 ). In this paper we report 
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on a similar phenomenon but we would like to point out that in principle it can be 
realized under other physical conditions. This possibility is based on the fact that a 
magnetic flux with time-dependence like in Eq. (|16[) may be realized with the help of a 
superconducting current. Indeed, according to resistively shunted junction model [20| 
the time dependence of the current going through the junction, which may produce 
the desirable external flux, reads 

Is(t)=I- j^'^-, , Z = ^R^~P C (20) 

I + l c coswi n 

where R is a resistor parallel to the junction, / is an externally applied constant circuit 
current and I c is a critical value of the persistent current. The time interval, during 
which the qubit is in state, depends on frequency u> from Eq. (116|) which according 
to Eq. (|20|) may be continuously tuned with an externally applied current /. 
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